
Goddard Space Flight Center 

GSFC STI Public Disclosure 
Export Control Checklist 


The Export Control Office requests your assistance in assuring that your proposed disclosure of NASA 
Space Administration '(NASA) Scientific and Technical Information (ST ) comply w„h 
the Export Administration Regulations (EAR, 1 5 CFR 730-774) and the International Traffic In Arms 
Regulations (1TAR, 22 CFR 120-130). The NASA Export Control Prograte .requires that 

and international presentation/pubhcation of Goddard Space Ce " ' [ ( ’ Scientific and 

through the GSFC Export Control Office in accordance with the NASA Foim 1676 NASA Sclent 
Technical Document Availability Authorization (DAA) process. Release of NASA information into 
nublic forum may provide countries with interests adverse to the United Stales with access to NAS 
technology Failure to comply with the ITAR regulations and/or the Commerce Department regulations 
may subject you to fines of up to $ I million and/or up to ten years imprisonment per violation. 

Completion of this checklist should minimize delays in approving most requests tor 

presentation/publication of NASA STI. 

Generally, the export of information pertaining to the design, development, production, manufacture, 
assembly operation, repair, testing, maintenance or modification of defense articles, i.e space g 
hardware ground tracking systems, launch vehicles to include sounding rockets and meteorological 
rockets, radiation hardened hardware and associated hardware and engineering units fo ^ hes ^^ e ^ e can 
controlled by the State Department under the ITAR. A complete listing of items covered by the TAR 
be accessed at http://lmd.gslc.na s a - g ov/exDort/home.htm- The export of information with «sP ec 
ground based sensors, detectors, high-speed computers, and national security and missile techno ogy 
items are controlled by the U.S. Commerce Department under the EAR- ^ the m formation mtended 
release falls within the above categories but fits into one or more of the following exemp 1 , 

information may be released. 


EXEMPTION I 

If vour information is already in the public domain in its entirety through a non-NASA medium and/or 
through NASA release previously approved by the Export Control Office the information is ; exempt from 
further review. If the information falls into this category, you may attest that you are using 
exemption by signing below. 


Signature 


Date 


EXEMPTION II 

If your information pertains exclusively to the release of scientific data, i.e. data pertaining to st udies 
of clouds soil vegetation, oceans, and planets, without the disclosure of information pertaining to articles 
controlled by the ITAR or EAR, such as flight instruments, high speed computers, or launch vehicles the 
information is exempt from further review. If the information falls into this category, you may attest 
you are using this exemption by signing below. 


Signature 


7 < 



Date 




i o «- 


GSFC 25-49 (1 1/01) Previous editions obsolete 




JOURNAL OF GEOPHYSICAL RESEARCH, VOL. 107, NO. 0, 10.I029/200IJA000507, 2002 


Long-term trends in interplanetary magnetic field strength 
and solar wind structure during the twentieth century 

I. G. Richardson 1 

NASA Goddard Space Flight Center, Greenbelt, Maryland, USA 


E. W. Cliver 

Air Force Research Laboratory, Hanscom Air Force Base, Bedford, Massachusetts, USA 

H. V. Cane 2 

NASA Goddard Space Flight Center, Greenbelt, Maryland, USA 

Received 28 September 2001; revised 20 November 2001; accepted 21 November 2001; published XX Month 2002. 

[il Lockwood et al. [1999] have recently reported a ~ 40% increase in t *\ e radial 
component of the interplanetary magnetic field (IMF) at Earth between 1964 and 1996. 
We argue that this increase does not constitute a secular trend but >s largely the 
consequence of lower than average fields during solar cycle 20 [ %, 6 ' m , 

comparison with surrounding cycles. For times after 1976 the average IMF strength has 
actually decreased slightly. Examination of the cosmic ray intensity, an indirect measure of 
the IMF strength, over the last five solar cycles (19-23) also indicates that cycle averages 
of the IMF strength have been relatively constant since ~ 1954. We also consider the 
origin of the well-documented increase in the geomagnetic aa index that occurred 
primarily during the first half of the twentieth century. We surmise that the coronal mass 
ejection (CME) rate for recent solar cycles was approximately twice as high as that tor 
solar cycles 100 years ago. However, this change in the CME rate and the accompanying 
increase in 27-day recurrent storm activity reported by others are unable to account 
completely for the increase in aa. Rather, the CMEs and recurrent high-speed streams at 
the beginning of the twentieth century must have been embedded in a background of slow 
solar wind that was less geoeffective (having, for example, lower IMF strength an or 
flow speed) than its modem counterpart. INDEX TERMS : 1650 Global Change; Solar vanab.Iity; 

2104 Interplanetary Physics: Cosmic rays; 2134 Interplanetary Physics: Interplanetary magnetic fields; 216. 
Interplanetary Physics: Solar cycle variations (7536); 7524 Solar Physics, Astrophys.es, and Astronomy. 
Magnetic fields; KEYWORDS: interplanetary magnetic field, geomagnetic activity, cosmic rays, so ar win , 
long-term trends 


1. Introduction 

[ 2 ] Recently, Lockwood et al. [1999] [see also Stamper et 
al 1999; Lockwood and Foster , 2000; Lockwood , 2001] 
reported a ~ 40% increase in the radial component of the 
interplanetary magnetic field (IMF) since the beginning of 
in-situ observations by spacecraft in the early 1960s. More- 
over, they used the geomagnetic aa index \Mayaud> 1972] 
and solar wind data to infer that a doubling of the solar open 
magnetic flux has taken place since ~ 1900. Here, we first 
examine IMF and galactic cosmic ray (GCR) data to 
determine in what sense, if any, the reported — 40% IMF 
increase from 1964-1996 constitutes a trend. The cosmic 
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ray data are relevant to this question because, as has been 
noted by several authors [e.g., Slavin and Smith, 1983; 
Cane et at., 1999a, 1999b; Belov, 2000], the GCR intensity 
tends to anticorrelate with the IMF strength. If there were 
a long-term increase in the IMF, then we might also expect a 
general decline in the GCR intensity with time. Such a 
decline has been reported by Stamper et al [1999] for the 
1964-1996 interval. In this paper, we suggest that these 
reported changes in the IMF and GCR intensity in 1964- 
1 996 are not indicative of true long-term trends. In partic- 
ular, IMF observations since 1976 actually suggest a decline 
in the average field strength. 

[ 3 ] As the second aspect of this study, we examine the 
cause of the remarkable rise in the geomagnetic aa index 
during the twentieth century, as noted by Lockwood et at 
[1999] and others [e.g., Russell, 1975; Feynman and 
Crooker, 1978]. The principal increase in aa was completed 
by ~ 1950 and average activity levels since then have 
remained relatively constant. The increase was so pro- 
nounced that the lowest activity levels observed after 
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1950 tend to be higher than the maximum levels observed 
early in the century. Since the long-term variability in aa 
closely reflects the envelope of sunspot maxima [e.g., Cliver 
et al., 1998] and the rate of coronal mass ejections (CMEs) 
tracks the sunspot number, at least in recent cycles [Webb 
and Howard, 1994; Cliver et al. , 1994], one possible cause 
of the secular increase in aa from ~ 1900 - 1950 is a 
change in solar wind structure over this period, with CMEs, 
which generate the vast majority of large geomagnetic 
storms [Gosling et al ., 1991; Richardson et al., 2001], 
occurring with increasing frequency. But was this the sole 
cause of the long-term increase in aa, or was there an 
increase in the geoeffectiveness of all components of the 
solar wind (i.e., slow solar wind, high-speed streams and 
CME-related structures)? We will address these questions 
by apportioning geomagnetic activity at the beginning of the 
twentieth century to the various solar wind flow types based 
on an analysis for recent cycles [Richardson et al., 2000, 
2002], and also by comparing geomagnetic activity for the 
aa minimum years of 1901 and 1977. 

2. IMF Variation During the Space Age 
2.1. Is There a Long-Term Trend in In Situ 
Observations of the IMF? 

[ 4 ] We first consider the report by Lockw’ood et al. [1999] 
and Stamper et al. [1999] that the radial component of the 
near-Earth IMF has increased in strength by ~ 40% since 
in-situ observations began in 1963. Figure lb shows Car- 
rington-rotation averages of the magnitude of the radial 
component of the IMF (\B r \) from the NSSDC OMNI data 
base. The solid line through these data shows the linear fit to 
annual IMF averages for 1964 1996 obtained by Stamper 
et al. [1999] (from their Table 2). Over this interval, the fit 
clearly indicates an increase in | B r \. We are reluctant, 
however, to interpret this increase in terms of a long-term 
trend in the IMF. Figure la shows Carrington-rotation 
averages of the solar 10 cm flux, an indicator of solar 
activity levels, for the three cycles (20-22) considered by 
Lockwood et al. and Stamper et al. together with the decline 
of cycle 18, cycle 19, and the rise of the current cycle (23). 
The radial component of the IMF clearly shows variations 
associated with the solar activity cycle which are most 
prominent in cycles 2 1 and 22. When the radio and IMF 
data are considered together, it appears that the rising IMF 
“trend” from 1964-1996 is primarily due to low field 
values measured during weak cycle 20, followed in cycles 
2 1 and 22 by a return to the higher IMF values inferred for 
cycle 19 from the intense 10 cm emission during that cycle. 

If fits are made to other intervals, then quite different trends 
are obtained. For example, the dotted line with a solid circle 
at each end in Figure lb shows that the fit for 1976-1996 
(i.e., omitting the first cycle of the Stamper et al. interval) 
has no overall trend. The fields observed from 1996 to 2000 
in cycle 23 appear to fall below the trend line obtained by 
Stamper et al. [1999] (though we note that the pattern in 
previous cycles (Figure 1) suggests that stronger fields 
might occur following the maximum of this cycle). In fact, 
a linear fit to data from 1976-present (the dashed line in 
Figure lb) indicates a decline in the average IMF strength 
during this period. Finally, a fit to the complete period of 
\B r \ data in Figure 1 (not included in the figure) shows an 


increase from ~ 3.03 nT to ~ 3.53 nT. This increase of 
~ 17% is much smaller than the ~ 40% reported by Stamper 
et al. [1999] for a subset of these data. Again it should not be 
interpreted as a “long-term” trend because it is sensitive to 
the endpoints of the relatively short period of data fitted. 

[5] The increase in the aa index (also ~ 40%) from 
1964-1996 reported by Stamper et al. [1999] is similarly 
sensitive to the interval considered. The aa index is shown 
in Figure 1c. As was the case for \B r \, the putative trend 
during the space age (solid line) appears to result from the 
weaker than (recent)average cycle 20. A fit to data for 
1954- 1996 (dashed line) indicates a much weaker upward 
trend of 0.09 ± 0.02 nT/year, while that for the complete 
interval shown in Figure 1 (not included in the figure) 
shows no significant trend (0.02 ± 0.02 nT/year). 

2.2. IMF Strength Since 1951: Implications From 
Galactic Cosmic Ray Observations 

[a] Another way of inferring whether there has been a 
change in the IMF intensity at Earth during recent solar 
cycles is to consider changes in the cosmic ray intensity. 
Figure Id shows solar-rotation averages of the cosmic ray 
intensity measured by the University of Chicago Climax 
neutron monitor (cut-off rigidity ~ 3 GV) since the begin- 
ning of observations in 1951. Note that the cosmic ray 
intensity is anti-correlated with solar activity levels, and that 
the size of the depression at solar maximum is roughly 
correlated with the size of the solar cycle and with the 
associated magnetic field increase, as has been previously 
reported [e.g., Slavin and Smith , 1983; Cane et al. , 1999a, 
1999b; Belov, 2000]. Stamper et al. [1999] noted a down- 
ward trend in the Climax counting rate from 1 964 to 1 996 
which they interpreted as evidence supporting a long-term 
trend in the IMF intensity. The solid line in Figure Id fitting 
the data for this period has a downward slope of ~ 0. 1 5 ± 
0.03%/year. This behavior results from the weak cosmic ray 
depression in cycle 20 followed by stronger depressions in 
the two subsequent solar cycles. Thus, as was the case for 
\B r \ and aa, we do not take the downward slope in the 
cosmic ray intensity during this three-cycle period as 
evidence for a long-term trend. If we consider instead the 
period from 1954 to 1996 (the minima preceding cycles 19 
and 23 respectively), the dashed line indicates that there is a 
negligible downward trend (-0.009 ± 0.02%/year). 
Because of the large amplitude of the solar cycle modu- 
lations compared to any long-term change over this interval, 
we note that the calculated long-term trend is extremely 
sensitive to the start and end times chosen for the fit. 

[7] Another way of examining the cosmic ray data for 
evidence of long-term trends is to compare the intensities 
during solar minima. The correct way to do so is to compare 
intensities in alternate minima where the solar global mag- 
netic field (A) has the same direction. This is because (see 
Figure 1) the intensity time profiles (peaked or flat when A< 

0 or > 0 respectively) and maximum intensity attained (at 
neutron monitor energies, higher when A< 0) are influenced 
by the differences in cosmic ray drift patterns in the helio- 
sphere which depend on A [e.g., Jokipii and Thomas , 1981 ; 
Kota and Jokipii 1983]. The counting rate of the data in 
Figure 1 has recovered to approximately the same level (to 
within < 20 counts/s, i.e. < 0.5 %) in minima of the same A 
epoch since 1954. This is a small fraction (~ 3%) of the solar 
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Figure 1 Comparison of (a) the solar 1 0 cm flux (a measure of solar activity levels, (b) the magnitude 
of the radial intetplanetary magnetic field component near the Earth, (c) the 

c i riQQQi Kiit this trend is strongly influenced by the weak cosmic ray 

reported by Stamper et a. [ L l954 , 9% indicates a negligible trend 

intensity dunng (alternate) solar minima (with the 
same ^direction of the solar global magnetic field A) shows no apparent trend. 
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cycle modulation, suggesting that if it is attributable to a 
change in the IMF intensity, this change is small compared to 
the solar cycle field variation. Thus, the cosmic ray data 
presented here show little evidence of a significant long-term 
trend (in particular a decrease) which might be consistent 
with a long-term change in the IMF strength. 

D] Before leaving this topic, we note that neutron mon- 
itor data counting rates may be influenced by factors such as 
changes in instrumentation (e.g. improvements/malfunc- 
tions) requiring careful renormalization, changes in station 
location, and long-term drifts in geomagnetic cut-offs 
because of temporal changes in the terrestrial magnetic 
field, the importance of which will depend on station 
location [e.g., Popielawska and Simpson , 1991; Shea and 
Smart, 2000, 2001], Thus, studies of long-term variations in 
cosmic ray intensity must ensure that the appropriate 
corrections have been made. In particular, a drift is known 
to have occurred in the Climax counting rate relative to 
other more stable neutron monitors and it is the renormal- 
ized data (J. A. Simpson and R. Pyle, “Post-1980 correc- 
tions to the Climax, Colorado neutron monitor counting 

rate”, http://ulysses.uchicago.edu/NeutronMonitor/Misc/ 
neutron2.htm!) that are used in Figure 1 . Ground level 
events (rare solar particle events which reach sufficiently 
high energies to be detectable by neutron monitors) have 
also been removed from these data. Nevertheless, these 
corrections are of minor importance to our main point here 
that the change in cosmic ray intensity inferred by Stamper 
et ai [1999] is the result of fitting to data from the weak 
cycle 20 and the subsequent two stronger cycles, and does 
not provide support for their claim of a long-term increase 
in the IMF during the same period. 


3. Increase in the Geomagnetic aa Index From 
~ 1900 to ~ 1950 

3.1. Was the aa Increase Due Solely to a Change in 
Solar Wind Structure? 

[9] We now turn to the second main topic of this paper - 
the cause of the secular increase in aa during the first half of 
the twentieth century. Figure 2b [after Vennerstroem, 2000] 
shows that, beginning in 1901, cycle-averages of the aa 
index exhibited a remarkable systematic increase during the 
first half of the twentieth century before settling at a higher 
level after ~ 1950. The aa index is correlated with rB s , 
where V is the solar wind speed, and B s is the southward 
magnetic field component [e.g., Crooker et al. , 1977; 
Crookerand Gringauz , 1993, and references therein’]. Thus! 
the increase in aa during the first half of the twentieth 
century could be caused by an increase in the average solar 
wind speed and/or B (assuming an approximately constant 
proportionality between B and B s ). 

[ 10 ] Alternatively, it might result from a change in the 
structure of the solar wind near the ecliptic plane, with, for 
example, the fraction of the solar wind populated by more 
geoeffective structures such as shocks and interplanetary 
coronal mass ejections increasing from ~ 1901 - 1 950 at 
the expense of less geoeffective high- and low-speed 
corotating streams. (Other possibilities, such as instrumen- 
tal or magnetospheric effects, seem to be insufficient to 
account for this increase in aa [e.g., Clilverd et al., 1998].) 
We will begin by considering whether the secular increase 
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Figure 2. Variation in the yearly averages of aa from 1868 
to 1999 (b), and for the 50 quietest days (a) and 50 most 
disturbed days (c) in each year [after Vennerstroem, 2000], 


in aa resulted entirely from a change in the solar wind 
structure. 

[n] First, we examine whether an increase in the rate of 
coronal mass ejections could have caused the long-term 
increase in aa . To do this, we use the results of Webb and 
Howard [1994], who combined coronagraph/photometer 
data from the Helios 1/2, Skylab, Solwind, and the Solar 
Maximum Mission spacecraft to obtain daily CME rates 
from 1974 to 1989, including the maximum of cycle 21. 
They showed that the CME rate tracked the solar cycle, and 
produced a correlation plot between yearly averages of the 
daily CME rate and the sunspot number (SSN). Assuming 
that this relationship holds for cycles covering a range of 
peak sunspot numbers, we can use it to infer the daily CME 
rate at the peak of cycle 14 (1902-1913). The smoothed 
peak (average) sunspot number in cycle 14 was 64 2 (32 3) 
compared with 164.5 (81.2) for cycle 21. From the relation- 
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Figure 3. D.stnbut.ons of the 3-hour ua ,ndex for (a) all solar wind regions ir ,1901 and^b) >for slow 
wind and (c) all solar wind in 1977. Arrows indicate average values, and the 1977 slow solar wina 
distribution is normalized to 1 year. Note that activity levels in 1901 were even less than those observed 
in slow solar wind in 1977. 


ship of Webb and Howard [1994], we then deduce a peak 
CME rate of ~ 1.1 CME/day for cycle 14. This is 50-60% 
lower than the peak - 2.5 CME/day rate for cycle 2 1 . Based 
on the analysis of Richardson et al [2000] [see also 
Richardson et al. , 2002], we determined that during the 
peak years of cycle 21, the Earth spent about one third of 
the time in each of the three principal components of the 
solar wind: CME-related flows (including transients, 
shocks, and post-shock flows), high-speed streams, and 
slow solar wind. Our CME-rate estimate for cycle 14 
suggests that Earth would have been embedded in CME- 
related flows for only ~ 15% of the time during the 
maximum of that cycle (assuming that these flows had 
similar characteristic sizes in both cycles). 


[ 12 ] Richardson et al. [2000] also determined the mean 
values of aa associated with these different types of solar 
wind flows during cycle 21. Taking average values (for the 
maximum of cycle 21; 1978-1982) of aa ~ 35 nT for CME- 
related flows, - 25 nT for high-speed streams, and - 15 nT 
for slow solar wind, a one-third mix of each solar wind type 
implies aa ~ 25 nT. If we reduce the fraction of CME-related 
flows to ~ 15% (appropriate for cycle 14) and assume that 
the remainder of the solar wind is equally divided between 
slow solar wind and high-speed streams, then we get an 
average aa of - 22 nT. Taking a more extreme case where 
the CMEs are embedded in slow solar wind flow yields aa ~ 
18 nT. However, Figure 2b indicates that aa early in the 
twentieth century was lower than each of these estimates (for 
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Figure 4. Comparison of the daily-averaged aa index in 1901 (shaded) and 1977, showing the 
persistently lower values attained in 1901. Some intervals of relatively enhanced activity in 1901 most 
clearly after day 200 and indicated by arrowheads, show evidence of recurrence at the solar rotation 
period, suggesting that corotating streams were present. 


example, the average aa during the maximum of cycle 14 
(1904 - 1909) was ~ 15 nT), showing that the increase in aa 
since this time did not arise solely from the inferred change 
in the CME rate. Only by completely excluding CMEs and 
assuming that only slow solar wind was present in the solar 
wind 100 years ago can we obtain a value of aa (~ 15 nT) 
matching that for the peak of cycle 14. However, we know 
that even great transient storms occurred early in the twen- 
tieth century (e.g., the September 1 909 storm associated with 
low-latitude aurorae [Silverman „ 1995]), so a CME-associ- 
ated component of the solar wind was certainly not com- 
pletely absent at this time. In addition, as we will show 
below, corotating streams were also present during this 
period. Thus, we conclude that the cause of the increase in 
aa from ~ 1900 - 1950 cannot be simply due to a change in 
solar wind structure. 

3.2. A Comparison of the aa Minimum Years of 1901 
and 1977 

[n] To examine the structure of the “background” solar 
wind at the beginning of the twentieth century, it is 
instructive to compare aa during a minimum year at the 
beginning of the twentieth century (1901) and for a recent 
cycle (1977). Distributions of 3-hour values of aa for these 
years are shown in Figures 3a and 3c, respectively, while 
Figure 3b illustrates the similar distribution for periods of 
slow solar wind in 1977 (normalized to 1 year). We note 
that the average aa for all solar wind in 1901 (~ 6 nT; 
shown by the arrow) was lower than that for slow solar 
wind in 1977 (~ 13 nT), which in turn was below the 1977 
value for all solar wind (~ 20 nT). Thus, the low levels of 
geomagnetic activity at the beginning of the twentieth 
century did not simply arise because the Earth was 
immersed at that time in solar wind with characteristics 
similar to those of slow solar wind detected during the space 
era. Eliminating both CMEs and high speed streams from 
the 1977 solar wind still leaves a solar wind that is twice as 
geoeffective as the average solar wind for 1901. Evidence 


that slow solar wind during 1977 is more geoeffective than 
that in 1901 can also be inferred from Figure 2 [after 
Vennerstroem [2000] where it is shown that the increase 
in average aa during the first half of the twentieth century 
(Figure 2b) was apparent on the quietest days (Figure 2c) 
presumably dominated by the slow solar wind. The activity 
increase was also seen on the most active days (Figure 2a) 
suggesting that CME-related structures likewise became 
more geoeffective during this period. 

[i4] Figure 4 compares time series of daily-averaged aa 
values in 1901 and 1977. Again, it is evident that activity 
levels on the quietest days in 1977 were higher than those 
generally found in 1901. Note that there were relatively 
brief intervals of enhanced activity in 1901. Some of these 
appear to recur at the ~ 27 day solar rotation period 
(examples are indicated by arrowheads), suggesting that 
they are associated with corotating high-speed streams [e.g., 
Crooker and Cliver , 1994]. The presence of recurrent 
structures is confirmed by the Lomb periodgram [Lomb, 
1976] for the 3-hr aa index in 1901 shown in Figure 5. 
Prominent peaks are centered on periods of ~ 26 - 28 days. 
These results show that the low values of aa at the 
beginning of the twentieth century do not correspond to a 
structureless ground state” of the solar wind. Figure 4 also 
shows that the recurrent activity increases in 1901 were less 
intense than those in 1977, indicating that the relative 
geoeffectiveness of corotating high-speed streams was also 
less in 1901 than in 1977. In fact, Clilverd et al. [1998] 
attribute the increase in aa during the first half of the 
twentieth century primarily to an increase in recurrent storm 
activity. 

[is] The reduced geoeffectiveness of slow- and high- 
speed streams ~ 100 years ago implies either a lower 
IMF strength at that time, or lower solar wind speeds, or 
both [Feynman and Crooker , 1978]. Since we obsen/e that, 
in recent solar cycles, variations of aa have been accom- 
panied by similar variations in the IMF strength (e.g., 
Figure 1), it is likely that the general increase in aa from 
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Figure 5. Lomb penodgram for the 3-hr averaged aa ,ndex in 1901 showing peaks in the power at 
approximately the whole and half solar rotation period (~ 27 and 13.5 days) indicating the presence o 
corotating structures. 


~ 1900 - 1950 involved a long-term change in the IMF 
such as has been proposed by Lockwood et al. [1999]. 
Moreover, since the solar-cycle variations in the IMF 
strength are manifested in each solar wind component 
[Richardson et al , 2000, 2002], it seems likely that a 
long-term increase in the IMF will involve high-speed 
streams and CMEs as well as slow solar wind. Note that 
because the long-term increase in aa occurred during the 
first half of the twentieth century, and aa has subsequently 
on average been more constant (Figure 1), a long-term 
increase in the IMF strength is not inconsistent with the 
absence of an increase during the second half of the century 
as inferred in Section 2. The large increase (nearly a 
doubling) of aa during the first half of the twentieth century 
indicates that changes in the IMF and/or solar wind speed 
were substantial during this period. 

4. Summary and Discussion 

[i 6 ] The Lockwood et al [1999] report that the Sun’s 
open magnetic field has doubled during the past century has 
generated a great deal of interest, leading to this special 
section of the Journal of Geophysical Research. Arge et al 
(this volume) point out that the related finding by Stamper 
et al [1999] of a 40% increase in the radial component of 
the IMF during the space age is not reflected in an increase 
in the open flux at the Sun deduced from measurements at 
three separate solar observatories since 1974. We show 
above that the inferred ~ 40% increase arises largely 
because of the difference in IMF strength between cycles 
20 (1964-1976) and 21 (1976-1986) (Figure 1), and 
should therefore not be regarded as a long-term or secular 


trend. For the post-1974 period considered by Arge et al, 
we find no evidence for an increase in the average radial 
component of the IMF. Indeed, there is a slight downward 
trend when one includes data for current cycle 23. Similarly, 
the ~ 40% increase in the aa index during 1964-1996 
reported by Stamper et al [1999] is largely the result of 
weak geomagnetic activity during cycle 20 and is not 
evidence for a long-term trend. 

[ 1 7 ] A consideration of galactic cosmic ray data, an 
indirect indicator of IMF strength, reinforces the evidence 
that there is no long-term increase in IMF since the early 
1950s (Figure 1). We note that Mursula et al [2001] 
inferred only a weak downward gradient in the average 
Climax neutron monitor intensity since 1951. They con- 
cluded that the trend found by Stamper et al [1999] and 
also by Lockwood [2001] and Lockwood and Foster [2000] 
combined with these authors’ assumption of a linear rela- 
tionship between the solar open magnetic field and the 
cosmic ray intensity, implies unphysically high cosmic ray 
intensities early in the twentieth century. Stozhkov et al. 
[ 2000 ] have also concluded that there is a weak long-term 
decline (few hundredths of a percent/year) in the cosmic ray 
intensity at a wide range of rigidities based on a more 
comprehensive study than we have attempted here. They 
suggest that this decrease could be explained bya supernova 
explosion at a distance of 30-150 pc about 10 —5x10 
years ago, rather than by changes at the Sun, interplanetary 
space, or the Earth’s magnetosphere, though they note that it 
is difficult to rule out such effects completely [cf. Ahluwalia 
and Lopate, 2001]. 

[is] As the second topic of this paper, we investigated the 
origin of the upward trend in aa during the first half of the 
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twentieth century. While the case for a clear trend in solar, 
solar wind, and geomagnetic conditions during the space age 
is problematic (Figure 1 ), the marked increase in aa between 
1900 and ~ 1950 (substantiated by corresponding increases 
in several other geomagnetic indices [Cliver and Ling, 2002] 
and its reflection in the sunspot number envelope [e.g., 
Cliver et al. , 1998]) provides strong evidence for a signifi- 
cant secular change in the Sun’s open magnetic flux as 
suggested by Lockwood et al. [ 1 999] and/or solar wind speed 
(see Feynman and Crooker [1978]). In this study, we enter- 
tained the idea that the low values o (aa for sunspot cycles at 
the turn of the twentieth century might have resulted from a 
reduction in the rate of CMEs and a circumstance where 
Earth experienced only slow solar wind, the least geoeffec- 
tive of the three basic flow types. We inferred that CMEs 
were present 100 years ago, but would have occurred at 
about half the present rate. This change in the CME rate, and 
an accompanying change in the frequency and extent of 
high-speed streams [Clilverd et al. , 1998] however are 
insufficient to account for the change in aa. 

[ 19 ] We presented evidence that the slow solar wind early 
in the last century was significantly less geoeffective than 
that of today [see also Vennerstroem, 2000]. We conclude 
that the rise in aa during the first half of the twentieth 
century did not result solely from a change in solar wind 
structure (e.g., an increased rate of CMEs) but also involved 
a general increase in the geoeffectiveness of the solar wind 
(related to increases in the IMF and/or solar wind speed) 
that is evident in the slow solar wind and most likely 
involved CMEs and corotating high-speed streams as well. 
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